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Abstract 

Exact results for the dynamic dimer and trimer structure factors of the one-dimensional s = 5 
XX model in a transverse magnetic field (|| z) are presented and discussed in relation to known exact 
results for the dynamic spin structure factors. In the framework of the Jordan-Wigner representation, 
the accessible spectrum of the dimer fluctuation operator is limited to two-fermion excitations whereas 
that of the trimer fluctuation operator involves two-fermion and four-fermion excitations. The spectral 
boundaries, soft modes, and singularity structure of the four-fermion excitation continuum as probed 
by the dynamic trimer structure factor are examined and compared with corresponding properties of 
the two-fermion excitation continuum, as probed by the dynamic dimer and transverse spin structure 
factors. 
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1 Introduction 



The theoretical and computational study of frequency-resolved quantum fluctuations and thermal fluctua- 
tions in many-body model systems is an important area of research for several reasons. Such fluctuations 
are observable directly or indirectly by a host of measuring techniques used in condensed matter physics 
and materials science. The shape of the spectrum (dispersions, bandwidths, gaps, soft modes, etc.), the 
spectral-weight distributions, and the singularity structure of the dynamic structure factors as measured or 
calculated for speciflc fluctuation operators yield detailed insights into the state of the material and reveal 
important clues about the susceptibility of the system to phase transitions with order parameters modeled 
after the fluctuation operators at hand 11. 

An increasing number of exactly solvable quantum many-body model systems turn out to be relevant 
for situations where certain degrees of freedom of a material are kinematically constrained to one spatial 
dimension. Many such situations offer the most detailed comparisons in many-body dynamics of exact 
theoretical results with direct experimental observation |21 El El ■ 

Various properties of dynamic spin structure factors of quantum spin chain models are observable in 
quasi-one-dimensional magnetic insulators, for example, via magnetic neutron scattering 5 . The properties 
of dynamic dimer and trimer structure factors, on the other hand, are important indicators of structural 
phase transitions driven by magnetic interactions such as in spin-Peierls compounds. Dimer fluctuations are 
key participants in phonon-assisted optical absorption processes of magnetic chain compounds and are thus 
observable in optical conductivity measurements |H1|71|S1|U1 ^1 . 

Here we consider the exactly solvable s = 5 XX chain with a magnetic field in the direction transverse 
to the spin coupling (in spin space) jlll I12| . The Hamiltonian reads 

N N 
n— 1 n— 1 

We set gfiB = l,h = 1, use the exchange constant as the energy unit, and measure the magnetic field Q in 
units of J. The Jordan- Wigner transformation to spinless lattice fermions maintains the bilinear operator 
structure, 

Af / 1\ 1 ^ 

n=l ^ ^ n=l 

A Fourier transform, = N^^/'^ ^^P (i'*") Cn, brings ()1.2|l into diagonal form: 

H = "^K^(clc^~]^ , A^^n + JcosK. (1.3) 

For periodic boundary conditions in Hl.l|l the allowed values of the fermion momenta depend on whether 
the number N j of fermions in the system is even or odd: K, G {(27r/iV) (n+ i)} \i N f IS even or S 
{(27r/A^)ri} if is odd. Fermion momenta within the first Brillouin zone are specified by integers n = 
-iV/2,-7V/2+ 1,. .. ,iV/2 - 1 (if iV is even) or n = -{N - l)/2,-(iV- l)/2 + 1, . . . , (iV - l)/2 (if N is 
odd). The Fermi level in the band A^ is controlled by the magnetic field fi. The number of fermions can 
vary between an empty band {N j = 0) and a full band (N f = iV) and is related to the quantum number S^' 
(the z-component of the total spin) of the same model in the spin representation: S^- = iV/ — A^/2. 

The dimer and trimer fluctuation operators for the XX model will be introduced in Sec. [3 The 
two-fermion and four-fermion dynamic structure factors associated with these fluctuation operators will be 



discussed in Secs.|21and^ respectively. Finally, in Sec. we give the conclusions and perspectives for future 
work. 



2 Fluctuation operators and dynamic structure factors 

Most fluctuation operators of interest are constructed from local operators An of the model system under 
consideration: 

1 ^ 

= — =^exp(iKn)v4„- (2.1) 



n=l 

Associated with each fluctuation operator (|2.1|l is a dynamic structure factor, 

SAA{n, u;)^27:J2 '^''^^'/^^'^ A |A) |' 6 {lu - Ex + Ey) , (2.2) 

AA' 

which describes fluctuations of a specific kind. Here E\ and |A) are the eigenvalues and the eigenvectors of 
H, and Z is the partition function. Of particular interest is the zero-temperature limit T = (i.e. /3 oo), 
where the thermal fluctuations fade away leaving pure quantum fluctuations in the wake. What remains in 
l|2.2|l are transitions between the ground state and all excited states that can be reached by the fluctuation 
operator A^: 

Saa{k,u) 2TrY,\{GS\A^\X)f S {lu ~ ux) , lox = Ex-Egs- (2.3) 

A 

The dynamically relevant spectrum observable in (|2.2|) or (|2.3|1 may vary considerably between fluctuation 
operators. Among other things, the spectrum is sensitive to their symmetry properties. 

For the s = ^ transverse XX chain the most important and most widely studied dynamic structure 
factors are those for the local spin operators 

Sn = cUn - ^, S+ = + isl = CXp ^171"^ c]cj ^ , s' = - IS^ = CXp ^-171"^ c]cj ^ C„. (2.4) 

At zero temperature the dynamic spin structure factor SzziK,i-u) is known to couple exclusively to the 
continuum of particle- hole excitations in the fermion representation, whereas Sxxii^,^) — Syy{K,uj) couples 
to excitations involving an arbitrarily high number of fermion excitations from the ground state |13l I14| . 

The fluctuation operators considered here are constructed from local spin operators on nearest and next- 
nearest neighbor sites. The dimer fluctuation operator _D„ and trimer fluctuation operator are obtained 
via (|2.1(l from 



and 



respectively. There is no unique way of defining dimer and trimer fiuctuation operators. The most suitable 
choice depends on the nature and the symmetry of the model system at hand. The operators 1)2.5(1 and 1(2. 6|l 
have the advantage that the associated dynamic structure factors Sdd{i^,oj) and Stt{i^,oj) can be analyzed 
exactly for the s — ^ transverse XX chain p. 1(1 in the fermion representation. 

As a motivation for the dimer and trimer operators used in this study, we offer a twofold argument. 
For a completely dimerized state, where nearest-neighbor spin correlations alternate between zero and a 



nonzero value along the chain, the operator ^/NDt^ plays the role of dimer order parameter. Likewise, for 
a completely trimerized state, where next-nearest neighbor spin correlations assume a period-three sequence 
of values zero, zero, nonzero, the operator ^fNT2Tr/z plays the role of trimer order parameter. 
Conversely, if we perturb the uniform XX Hamiltonian by interactions of the form 

N 

Hd =e^cos {■Kn)Dn (2.7) 

n=l 

or 

i^T cos (^y 71 jr„, (2.8) 

the ground state becomes dimerized or trimerized, respectively. In the former case, nearest neighbor cor- 
relations are modified by period-two perturbative corrections of order — e, -f e and in the latter case by 
period-three corrections of order — ^e, — jS, +£• 

We may formally introduce the polymer fluctuation operator of order /, "P^^ via H2.1|l from 

+ (2.9) 

It includes the dimer and trimer operators for I — 1,2, respectively: Vn^ = -D„, Vn^ = T„. From the fermion 
representation of the polymer operator H2.9|l as carried out explicitly in (|2.5() and (|2.6() for the lowest two 
orders, it is evident that the dynamic polymer structure factor S-p-p{k,lu) at zero temperature will involve 
2TO-fermion excitations with m = 1,2, ... ,1 from the ground state. For an infinitely long chain (N —>■ oo) 
the polymer fluctuation operator and the function S-p-p{k,lu) may thus serve useful roles in attempts to 
understand the enormously complex dynamic spin structure factors Sxx{n,u!) — Syy{K,uj). Such tools by 
which the complexity of the dynamically relevant excitation spectrum can be gradually and systematically 
increased are not only useful for the calculations in the fermion representation as performed here but also for 
the recently developed techniques of calculating transition rates for the XX model in the framework of the 
Bethe ansatz |15l 1161 ITT) . The time-dependent polymer correlation function is related, in the limit Z — s- oo, 
to time-dependent spin correlation functions as follows: 

{V^nKt)v'llM) ^'^^ 2«(i)<+„(0))2 + 2(.-(t)4+„(0))2. (2.10) 
Note that (s^(i)s^_,_^(0)) is nonzero only if ^ 0. 



3 Two-fermion dynamic structure factors 

We start with the dynamic quantities which are governed by particle-hole excitations. The equilibrium 
time-dependent correlation functions for the operators sfi(t) and D„{t) can be evaluated directly: 

« W<+i(0)) - {s^r = ^T. e-'^-^--^^' exp (i [A., - A.J t) n^, (1 - n.J , (3.1) 

Ki ,K,2 

{D^{t)D^+i{0)) - {Df = ^ E ^=0^' ^l±^e-i(-^-«^)' exp (i [A., - A.,] t) (1 - J , (3.2) 
where = 1/ (1 + exp {f3A^)) is the Fermi function and 

(^^) = ^E«) = -2^Et-h^' (3-3) 

n— 1 K 



1 ^ 1 RK 

= ]^E(^«) --^E^o^'^tanh^. (3.4) 

n— 1 K 

The associated dynamic structure factors, 

Saa{>^,oj) = ^exp(-iKO / dtexp(ic^i) ((A„(t) - {A)) (A„+KO) - (A))), (3.5) 
1=1 

all of which involve two-fermion transitions, are obtained by Fourier transform. The resulting expressions 
for — > oo can be brought into the form 



( dKin^i (1 - n«i+«) 5{uj + A^, - A^i+^) = V — TTT^TTT' 
j-TT 2| Jsm I cos (I + K*j I 



■r-^ cos^ (^ + k"^) n^* (1 — 1 K*) 

5z,z,(«, c) = ^ Y ■ ^ .cr ' (3-^) 

2| Jsm ^ cos (I + K*j I 
where — tt < k* < tt are the solutions of the equation 

= -2 J sin I sin + K*) . (3.8) 

The dynamic structure factors H3.6() . H3.7|l are governed by the two-fermion (particle- hole) excitation 
continuum the properties of which were examined in Refs. |18l This continuum is well visible in Figs. 
?? and ?? below, which we include mainly for comparison with new results. At zero temperature, T ~ 0, the 
two-fermion excitation continuum has the following lower, middle and upper boundaries in the (k, a;)-plane^ 



- — - = 2 sm — 
J 2 



— = 2sm-sml 



sin (f - «) I , (3.9) 
3in(^-^a), (3.10) 
2 sin — sin ( — ha), if < k < tt — 2a, 

Yl-^ I ' " " (3.11) 

"^1 2sin — , if TT — 2a < K < TT, 

\ 2 - - > 

respectively. The parameter a = arccos (f2/| J|) varies from tt when Vl = — |J| to when — \J\. The 
w-profiles at fixed k of the two-fermion dynamic structure factors may exhibit square-root divergences (a 
common density-of-states effect in one dimension) when oj 2\J\ sin(K/2). At T > the lower boundary of 
two-fermion excitation continuum is smeared out. The spectral weight in H3.6|l and (|3.7(l is now confined to 
\uj\ < 2|J|sin(K/2). 

Closed- form expressions for the two-fermion dynamic structure factors (|3.6|) . I|3.7|l exist in the low- 
temperature and high-temperature limits. At T = we have 



Q{uj — uji)Q{ll!^ — to), if < K < TT — 2a, 

^4 J2 sin^ f - [ {e{Lj-LJi) + e{Lj-Lj^))e{uJu-uj), if 7r-2a<K<7r, 



(3.12) 



/4J2sin2f-^2 j 0(^_^^)0(^^_^)^ if o<K<7r-2a, 

4j2sm''f {e{uJ-UJl)+e{uJ-LO,n))&iiOu-Uj), if TT - 2a < K < TT, 



^We assume that < k < tt in the rest equations of this Section; these equations are valid also for — tt < k < after the 
change k — » — k. 



and at T —^ oo we have 



u) - , ^ =e (21 J| sin 5 - H) , (3.14) 



2^/4J2 sin^f 



5z,i,(«,c.) = ^——-4^ e (2|J|sin- - Ic^l) . (3.15) 

The zero-tcnipcraturc rcsuhs for Szz{k,uj) can be found in Eq. (2.3) of Rcf. 18^ and for Sdd(k,uj) at = 
in Eq. (3.2) of Ref. 0. 

In Figs. ?? and ?? we show the dynamic structure factors (|3.6|l . (|3.7|) at zero temperature T — and 
different values of the transverse field (panels a, b, c), and at T — > 00 (panel d). The results for T 00 are 
independent of 17. As we can see, the two-fermion dynamic structure factors are nonzero within the two- 
fermion excitation continuum in the (k, CL;)-plane. Their spectral-weight distributions arc controlled by the 
Fermi functions, the multiplicity of the solution of Eq. I|3.8f) . the singularities in the density of one-particle 
states, and the explicit form of the rest of the integrand in (|3.()|) . H3.7|l . Another two-fermion dynamic 
quantity will be presented in Sec. 0] namely the two-fermion contribution to the dynamic trimer structure 
factor, s!^^{k,uj). 

4 Four-fermion dynamic structure factor 

Next we consider the dynamics of the trimer fluctuations. The method remains the same but its execution 
is more tedious. In addition to two-fermion transitions also four-fermion transitions contribute to the trimer 
fluctuations. The expression for the equilibrium time-dependent trimer-trimer correlation function reads 

(r„(i)r„+KO)) - = ^J2 ^^'^ ^^2) e~'^^^~^^^' exp (i [A«, - A«,] t) (1 - n.J 

ki ,k2 

+ 774 E C(4^('«i,«2,«3,«4)e~i('^^+''^-«^-''^"exp(i[A., -l-A,,-A,3-A,Ji) 

X Ui^^n^^ {1 - n,^^) {1 - n^J , (4.1) 

where 

1 ^ 

{T)^^Y.^Tr,)^C2 + 2cl-2coC2, (4.2) 



71=1 



C(2) (ki, K2) = (1 - 2cof cos^ (ki + K2) + 4ci (1 - 2co) (cos^ + y) + cos^ -|- K2)) 
+4:cl (cos^ Ki + cos^ K2) + 8 (-C2 + cl + 2C0C2) cos^ ^ + 8cJ cos^ ^ 



.,2 

2 

-t-4c2 - 8ci - 8c^ Acj + WcqCi - BcqCs + 16ciC2, (4.3) 



/ 2 '^i 

-4ci (1 — 2co — 4c2) ycos — + COS" 



^ (ki, K2, K3, K4) = 16sm sm cos . (4.4) 

Here we have introduced the function Cp — (l/N) cos (pn) n„. For N 00 and at zero temperature we 
have Co = 1 if < Co = a/n ii -\J\ < n < \ J\, co = if |J| < Cp = (-sgn( J))^ sin (pa) / (p7r) if 



\n\ < \ J\ and Cp = otherwise {p — 1,2, . . . ). At T ~> cxd we have Cp — ^(5p.o- The resulting dynamic trimer 
structure factor (|3.5|) ior N oo then has the following form: 

Stt{i<-,^) = s!^^{k,u!) + s!^^{k,uj), (4-5) 

where 



s'^^{k.,uj) = I dKiC*^' (ki,ki + (1 - Ti^i+h;) 5 (cj + A^i - Ak^+k) , (4.6) 



1 f 

SI^^{k,Uj) = ^— J J dKl J dK2 J dK^C'^^'' {ki,K2,K3,Ki + K2 - K3 + k) 

X (1 - n^^g) (1 - 11^^+,^^^^,^+^) S {u + A^^ + A^^ - A^,, - A^^+i^,^^f,^+^) . (4.7) 

The spectral weight in this quantity comes from both the two-fermion (one particle and one hole) excita- 
tion continuum and the four-fermion (two particles and two holes) excitation continuum. Let us first discuss 
the properties of the four-fermion excitation continuum and then the properties of Stt{i^t^)- At T = the 
four-fermion excitation continuum (for J = — \J\ < 0) is determined by the conditions 



■j — r — cos Ki -|- cos K2 — COS K3 — COS K4, K = —Hi — K2 + K3 -|- K4 (mod(27r)), 

n n n n , , 

COSKi>— , C0SK2>j-J|, C0SK3<— , COSK4<|-j|, (4.8) 

— TT < Ki_2,3 < TT, — TT < K < IT . Equatious H4.8|l imply that the four-fermion excitation continuum (like 
the two-fermion excitation continuum) exists only if the magnetic field does not exceed the saturation field: 

\n\ < \j\. 

An analytical expression for the lower boundary of the four-fermion excitation continuum in the (k, uj)- 
plane depends on Q and k and is given by one of the following expressions (see Appendix for additional 
details): 

^l^'' „ . |K| . / \K 



2sinYsin U-Y'' ^^'^^ 



(2) 



= 4cos-cos( a+ Y ) , (4.10) 



^ = -2sin(a + ^jsin(2a + ^), (4.11) 



^^-2sm(«-M)..(,_M,, (,,2) 



|J| 4 V 4 

The range in k over which for a given H. one of the expressions l|4.9|l - 14.13|l forms the lower boundary of 
the four-fermion continuum can be read off Fig. ??. The darkness in this gray-scale plot is a measure of 
the size of the energy threshold of the four-fermion continuum (white means zero excitation energy i.e. a 
soft mode) . The boundary between the region i (where cu^^^ is the lower boundary) and the region j (where 



w^'''' is the lower boundary) follows from the matching condition cj^^''* — ajp^ and is given by the formula 
|«;| = lij{a) where 



, / X , tan a — V tan^ a — 3 , , 2tt 

Ii2{a) = 4arctan , \k\ < (4.14) 

Zi3(a)=^-a, |<|K|<y, (4.15) 

/14(a) = 2a, (4.16) 

/23(a) -27r- 4a, (4.17) 
1 27r 

/34(a) = |k| +cosa- -, Y<|K|<7r, (4.18) 

/45(a) = 4a. (4.19) 
The boundary between regions 2 and 4 is determined by the cubic equation 

sin |k| — 2 sin -Li^ tan"^ a + ^3 + 2 cos + 3 cos tan^ a 

\k\ \ K 

2sin-y- + 3sin|K| j tana + 3 + 2 cos - - cos k 0. (4.20) 

Typical lower boundaries of the four-fermion continuum for several values of fl can be seen in Fig. ??. The 
soft modes according to (|4.9|) - 14.13|l are given by 

|ko| = {0,27r-4a,2a,4a}. (4.21) 

Alternatively, the soft modes 14.21() may be determined directly from H4.8|l . They occur when coski = 
cos K2 = cos K3 = cos K4 — cos a. 

The upper boundary of the four-fermion continuum for < < | J|/\/2 is given by 

^^4c«^, ,4,22) 

For \ J\/V2 <n<\J\ the upper boundary is given by H4.22|l only if |k| < 4a whereas, if 4a < |k| < tt, it is 
given by 

(2) / I 



|j| 4cos-cos(^a-Yj (4.23) 

(see Appendix) . The upper boundaries of the four-fermion continuum for several values of CI can be seen in 
Fig. ?? in comparison with the corresponding two-fermion continuum. 

The four-fermion continuum always contains the two-fermion continuum. The lower boundaries coincide 
in part. The upper boundaries are different. In the zero field case we have uji = | J| sin |k| for both continua. 
The upper boundaries are lUu ~ M-A cos j and w„ = 2| J| sin ^ for the two-fermion and four-fermion continua, 
respectively. As the saturation field f2 = | J| is approached from below, the two-fermion continuum narrows 
to a branch and then disappears whereas the four-fermion continuum remains an extended region, bounded 
by uji — 4| J| sin^ j and uju = 4| J| cos^ j, and then disappears abruptly. 



Now consider the equation 



. [ [COSKI + cos K2 — cos K3 — cos (k + Ki + K2 — K3)] ) =0. (4.24) 



It is satisfied for 



|J| -2sinY, (4.25) 



^=4sinM, (4.26) 



(3) 

UJ s K, , ^ 

— =4eos-. (4.27) 



Thus, for K or w approaching the curves (|4.25() - H4.27II in (k, a;)-space, the quantity 



S{k, bj) 



/TT />7r />7r 

d/ii / dK2 / d/iaS* (ki, K2, K3, k) 
-TT J —TT J —IT 

X (5 (cj — I J| COS Ki — I J| COS K2 + I J| COS K3 + I J| COS (k + Kl + K2 — K3)) (4.28) 



exhibits cusp singularities (akin to dcnsity-of-states effects in three dimensions). The exact nature of the 
cusps also depends on the factor 5 (ki, K2, K3, k), which varies between different dynamic structure factors 
with a four-fcrmion part. It always includes the factor 71^-^11-^2 (1 ~ "ks) (1 ^ '^Ki+Ka-Ks+K) as can be seen in 
expression (|4.7|l . In Fig. ?? we show the cj-dependence of S{k,ijj) as given by H4.28|l at k = 0, ^, tt 
when S {ki, K2, K3, k) — 1 and when S (ki, K2, K3, k) = ni^^Uf^^ (1 ~ '^kh) (1 ~ 11^1+^2-^3+^) fo'" several values 

of r2. 

At T > the lower boundary of the four-fermion excitation continuum is smeared out and the upper 
boundary becomes = 4| J| cos ^. 

The properties of the multimagnon continua of quantum spin chains have been examined in some detail 
in the recent paper of Barnes |20| . In particular, the lower/upper boundary of the two- and higher magnon 
continua were determined. It was shown that the boundary curves under certain conditions may exhibit 
discontinuous changes in composition and cusps. Moreover, a behavior of the density of (two- and higher 
magnon) states on the continuum boundaries and within the continuum was considered and the existence 
of discontinuities was pointed out. These features of one-dimensional quantum spin systems are expected 
to become accessible experimentally in high-resolution inelastic neutron scattering on alternating chain and 
ladder materials. 

Interestingly, the s — ^ transverse XX chain which can be mapped onto nonintcracting spinless fcrmions 
provides an excellent example of a system whose dynamic properties are governed by continua of multifcrmion 
excitations. In particular, the dynamics of trimcr fluctuations provides a direct motivation for analyzing the 
four-fermion excitation continuum. Unlike in the analysis reported in Rcf. 20 , where the statistics of the 
elementary excitations (magnons) is not known, here the quasiparticles are known to be fermions and the 
consequences are fully taken into account. 

Finally, let us examine the explicit expression for the dynamic trimer structure factor S'tt('*, w) l|4.5(l . 
In Figs. ??, ??, ??, ?? we present the zero-temperature dynamic trimer structure factor at different values 
of f2. In Fig. ?? we present the same quantity at infinite temperature. We show separately the two- 
fermion contribution (panels a) and the four-fcrmion contribution (panels b) as well as the sum of these 



contributions (panels c). We observe how the spectral weight is spread across the four-fermion continuum. 
We also see that the two-fermion contribution stands out in terms of spectral weight. The two-fermion and 
four-fermion contributions are comparable in intensity at T = and small fi. As increases the two-fermion 
contribution becomes more important and it completely dominates as 51 ^ |J|. In the high-temperature 
limit the two-fermion contribution is very dominant but the four-fermion continuum is still in evidence. 

5 Conclusions 

In summary, we have investigated some aspects of the dynamics of the s = transverse XX chain examining, 
in particular, the dynamics of dimer and trimer operators. For this purpose we have calculated several 
dynamic structure factors on a rigorous basis within the Jordan- Wigner representation. While the dynamic 
dimer structure factor S'dd(k,cj) and the dynamic transverse spin structure factor Szz{k,ll!) are governed 
by fermionic one-particle-one- hole excitations, the dynamic trimer structure factor STTiK,Lu) also contains 
contributions from two-particle- two- hole excitations. We have described the structure of the two-fermion 
and four-fermion excitation spectra in some detail and investigated the distribution of spectral weight in 
Stt{k,uj) across these continua at zero and nonzero temperature. In particular, we have established the 
boundaries of the four-fermion spectral range, the locations of soft modes, and the singularity structure, 
which includes one-dimensional and three-dimensional density-of-states features (van Hove singularities). 

An alternative technique to evaluate dynamic structure factor of quantum spin chains is based on the 
Bethe ansatz solutions I^IEIEI- Recently such an approach has been applied to the s — ^ XX chain 
Moreover, the relation between spinous or magnon quasiparticles and Jordan- Wigner fermions was discussed 
in some detail. It will be interesting to interpret the two-fermion and four-fermion excitations discussed here 
in terms of the Bethe ansatz solution as studied in Ref. 
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Appendix 

To find the lower and upper boundaries of the four-fermion excitation continuum (|4.8|l at fixed and 
—IT < K < t: we search for the extrema of cj/j J| as given in 1)4. 8|l and the values of ki, K2, ks, at which 
such extrema occur. Typical results are reported in Fig. ?? (lower boundary) and Fig. ?? (upper boundary). 

In Fig. ?? we show the dependence on |k| of ki, K2, K3, where J| assumes a local minimum. We 
distinguish five different regions. The global minimum yields the lower continuum boundary, li Q < k < Ka 



Ki = a — K, K2 — K3 — K4 — a 



(A.l) 



and 

(1) 

— COS (a — k) — COS a — 2 sm — sm [^^12) ~ TjT' i^-^J 



if Kq < K < 



and 



Ki = K2 = a, K3 = K4 = — + a — TT (A. 3) 



(2) 

CJ; f n \ K / K\ U), 

- — - = 2 cos a — 2 cos — h a — tt =4 cos — cos a H — = -rrr (A. 4) 
|J| V2 / 4 V 4/ |J| 

etc. The values of Ka, Kfc, Kc, 'td follow from the matching conditions. 

In Fig. ?? we show the dependence on |k| of ki, K2, K3, K4 where (j^/I J| assumes a local maximum. We 

distinguish two different regions. The global maximum yields the upper continuum boundary. If < < ha 

Kl = K2 = --, K3 = K4 = -7r+- (A. 5) 



and 



ii KA < K < TT 



and 



(1) 



Ki = K2 = —a, K3 = K4 = —a — TT + — (A. 7) 



(2) 

Wu„ n / , 1^ / '^u /A 

2 COS a + 2 cos a —4 cos — cos a — - — - ( A.i 

V 2/ 4 V 4/ J 



1^, 

etc. From the matching condition we find ka — 4q; 
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FIGURE CAPTIONS 



Fig. 1. S^^{k,uj) at T = and (a) O = 0, (b) f7 = 0.3, (c) il = 0.9, and (d) at T ^ oo (independent of 
il; only w > is shown). 

Fig. 2. Sdd{k,uj) at T = and (a) Q = 0, (h) n = 0.3, (c) Q = 0.9, and (d) at T ^ oo (independent of 
f2; only w > is shown). 

Fig. 3. Lower boundary uji = min (^^P''^ i J — I7 • ■ • 1 5 of the four-fermion excitation continuum versus 
wave number n and transverse field Q (for | J| = 1). 

Fig. 4. Lower and upper boundaries of the two-fermion and four-fermion continua for | J| = 1 and f2 = 0, 
fi ~ 0.3, fl — 0.6, fl — 0.9. The two-fermion continmmr is shown shaded. 

Fig. 5. S{k,lj) as given in l|4.28(l versus ui at k ~ 0, ^, tt with S {ki, K2, k^, k) = 1 (bold curves), 
and S {ki, K2, K3, k) = rii^-^n^^ (1 — n^g) (1 — n^^+Kg-Ks+K) for = (solid curves), fl = 0.3 (long-dashed 
curves), fl — 0.6 (short-dashed curves), = 0.9 (dotted curves). Vertical lines mark the values of wP'', 
j — 1, 2, 3 as given in H4.25|l - (|4.27(l . Note the different vertical scales left and right. 

Fig. 6. STTiK,uj) at T = and 17 = 0. Separate plots are shown for (a) s!^^{k,ll!), (b) s!^^{k,uj), and 
(c) the sum Stt{k,uj). 

Fig. 7. Stt{i^,^) at T = and — 0.3. Separate plots are shown for (a) ^^^^(k, w), (b) S^^{k,uj), and 
(c) the sum STrin^uj). 

Fig. 8. Stt{k,uj) at T = and = 0.6. Separate plots are shown for (a) s!^^{k,uj), (b) s!^^{k,lu), and 
(c) the sum Stt{>^7^)- 

Fig. 9. Stt{h,i^) a,t T — and = 0.9. Separate plots are shown for (a) S^I{k,uj), (b) s!^^{k,uj), and 
(c) the sum Stt{h,^)- 

Fig. 10. Stt{i^,^) at T — > 00 (independent of il; only w > is shown). Separate plots are shown for (a) 
s!^^{k,lu), (b) s!^^{k,uj), and (c) the sum STriK^Lu). 

Fig. 11. Search for the lower boundary of the four-fermion excitation continuum. Shown arc the values 
of Ki, K2, K3, at which a minimmn of tij/l J| as given in H4.8|l occurs at f2 = 0.3| J| and — tt < k < tt. The 
dependences ni and K3 on k are shown by dashed curves, the dependences K2 and K4 on k are shown by 
dotted curves, the dependence of the minimal value of tj/l J| on k is shown by solid curves. 

Fig. 12. Search for the upper boundary of the four-fermion excitation continuum. Shown are the values 
of Ki, K2, K3, K4 at which a maximum of J| as given in H4.8|l occurs a,t fl = 0.9| J| and — tt < k < tt. The 
dependences of ki(= K2) and K3{— K4) on k are shown by dashed curves, the dependence of the maximal 
value of w/l J| on k is shown by solid curves. 
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